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We consider a supersymmetric U(N) gauge theory with matter fields in the adjoint, funda- 
mental and anti-fundamental representations. As in the framework which was put forward 
by Dijkgraaf and Vafa, this theory can be described by a matrix model. We analyze this 
theory along the lines of [F. Cachazo, M. Douglas, N.S. and E. Witten, "Chiral Rings and 
Anomalies in Supersymmetric Gauge Theory" [hep-th / 02 1 1 1 70|| and show the equivalence 
of the gauge theory and the matrix model. In particular, the anomaly equations in the 
gauge theory is identified with the loop equations in the matrix model. 
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1. Introduction 

Recently Dijkgraaf and Vafa J1] were motivated by earlier work ||-^§ to conjecture an 
interesting relation between SUSY gauge theories and matrix models. Many authors have 
added matter in the fundamental representation ||-25[] to this framework. We will examine 



a general theory of this form which includes the various examples of 0-25] as special cases 



We will follow the point of view of pq] , and will extend it to this case with fundamental 
matter. 



We consider an M = 1 supersymmetric U(N) gauge theory with matter in the adjoint 
Vf fundamentals Q* and N 
The tree level superpotential is 



$, Nf fundamentals Q* and Nf anti-fundamentals Qj (/ and / are the flavor indices). 



W tree = Tr W($) + Qjm f f (<S>)Q f (1.1) 

where we suppressed the color indices. The function W and the matrix m are taken to be 
polynomials 

WW = E k^T^ i+1 

. " . < L2 > 

ml(z) = \^ ml ,z h 1 



fe=i 



We consider arbitrary N and Nf. For nNf > N these theories exhibit a certain 
duality f27j - |2~9[1 exchanging <-> nNf — N, and we hope that this discussion will shed light 



on it. 

In section 2 we analyze this quantum field theory focusing on its chiral ring, the 
anomaly equations and the low energy effective superpotential. In section 3 we consider 
the corresponding matrix model and prove its equivalence to the gauge theory. 

2. Gauge theory considerations 

2.1. Anomaly equations 

As in |?6| we will be interested in the chiral operators 



4* z-* , , 

1 WJJV a 
R( z ) = E_Tr a 

MUz) = Q 7 —^-rQ f 
f f z - $ 
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Since W a Q f and QjW a are not in the chiral ring, there is no need to include more opera- 
torJ. Also, since the gauge group is U (N) rather than SU (N) we do not include "baryonic 
operators." 

We perform the following five independent transformations in the functional integral 

5$ ' 



5$ 



z - $ 
z - $ 



S $ = V I^_ (2.2) 
z — $ 

SQf = j^X f r Q f ' 
5Q T =\ f ~Q Tr ' 



'/ f^f'z-Q 



where Aj-, and A— are constant flavor matrices. The anomaly equations of these transfor- 
mations are 

Tr 4f + Q 7 ^-_iQ/ = 2R(z)T(z) + w a (z)w a (z) 

z — J z — $ 

47T 2 - $ 

o^r ' =R(z) (2 1) 

32tt 2 v ; V' 6 ) 

where the equalities are only in the chiral ring; i.e. they are up to SUSY commutators. The 
0(1/ z) term in the first equation is the equation of motion of Tr $. The 0(1/ z 2 ) term 
in that equation and the 0(1/ z) terms in the fourth and fifth equations are the Konishi 
anomaly I^JtH] for Q and Q respectively. The other equations and the other powers of 



1/z are various generalizations. 



1 We thank E. Witten for a useful discussion on this point. 
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Since and m^($) are polynomials, we can replace these equations with 

[W'(z)T{z)]_ + tr [m'(z)M(z)]_ = 2R(z)T(z) + w a (z)w a (z) 

[W'(z)w a (z)]_=2R(z)w a (z) 

[W'{z)R{z)]_=R(zf 



(M(z)m(z)) f f ' 



R(z)5 f f ' 



(2.4) 



(m(z)M(z)) 



f 



R{z)5 



f 



Here m(z) and M(z) are matrices in flavor space. We multiply such matrices as (AB)^,, = 
Af,Bf„ and tr denotes a trace over the flavor indices. In the last two equations we used 
the fact that the corresponding equations in (|2.3| ) are satisfied for every A and A. These 
identities are also true for the expectation values of the operators (|2.1| ) . 
Writing the third equation in (|2.4j) as 



(2.5) 



W\z)R{z) + \f{z)=R{zf 
with a polynomial f(z) = —4 [W'(z)R(z)} + of degree n — 1, we solve for R(z) 



2R(z) = W'(z) - ^W'(z) 2 + f(z) 



(2.6) 



Similarly, we write the last two equations in (|2.4| ) as 

(M(z)m(z) - q(z)) f f ' = R(z)5 f f ' 
(m(z)M(z) -q(z)) f 7 = R(z)d f 7 



(2.7) 



in terms of polynomial matrices 

q(z) = [M(z)m(z)] + 

(2.8) 

q(z) = [m(z)M(z)} + 
We can now solve for the matrix M(z) 

M(z) = R(z)m- 1 {z) + q(z)m- 1 {z) = R(z)m- 1 {z) - m- 1 {z)q{z) (2.9) 

where m~ 1 (z) is the inverse matrix. For consistency the polynomial matrices q(z) and q(z) 
are related by q = mqm -1 . Using the solution for R(z) in ( |2.6| ) the expression for M(z) 
( pD| ) leads to a solution for M(z) in terms of the polynomials f(z) and q(z). 



The polynomial f(z) is kept arbitrary, but the polynomials q{z) are fixed as follows. 
From ( |2.9|) it is clear that M(z) is singular when m(z) has a zero eigenvalue. The freedom 
in the polynomials q(z) is exactly such as to remove these singularities. To see this, consider 
for simplicity the case Nf = 1, where the flavor matrices are one dimensional. m(z) is 
a polynomial of degree / ( |1.2|) , and therefore it has I zeroes zj. For large z the resolvent 
M(z) behaves as \ and therefore q(z) = [M(z)m(z)} + is a polynomial of degree / — 1. The 
I coefficients of q can be tuned to set the I residues of M(z) at zi to zero. 

Now that we have solved for R(z) and M(z), it is simple to solve for w ot (z) and T(z) 
using the first two equations in (|2.4|). In doing that we need to introduce new polynomials 



(2.10) 



Pa (z) = -4[W'(z)w a (z)] + 

c(z) = -4[W'(z)T(z)} + 



Note that the solution for R(z) is as in the theory without the fundamental matter, but 
T(z) is different; it depends on m(z). 

2.2. The effective superpotential 



We parametrize the superpotential of the adjoint field as 



W'{z)=g n \[{z-a i ) 



i=i 



in terms of its stationary points a^. We consider the classical vacuum 



(2.11) 



(<&> 



0'2 



a 3 



\ 



J 



(2.12) 



where ai occurs Ni times. We assume that VF($) and m($) are generic and in particular 
m(ai) does not have a zero eigenvalue. Since all the matter fields acquire masses, the low 



energy degrees of freedom are only the gauge fields in Y\ i U(Ni). Following [26] we express 
them as the gauge invariant objects 

Si = — <£ R(z)dz 

1 

2tH 

1 
2mi 



c, 



= 7T-. $ w a (z)dz (2.13) 



c, 



N, = — <b T(z)dz 



Ci 

where Cj is a contour around dj. 

Since we assumed that m is generic, all the quarks are massive and are integrated 
out. Therefore, our effective Lagrangian will not depend on the "meson operator in the % 
group" 

M (l) = — I M{z)dz (2.14) 
2ttz J Cl 

As in |2(J (see also |3"2]j ) we integrate out the massive fields in perturbation theory. 
It is straightforward to add the fundamental matter to the discussion of the Feynman 



diagrams in [p6 |. Assume for simplicity that the group U(N) is unbroken; i.e. N\ = N. 



Then considerations of symmetry and holomorphy along the lines of |33| constrain the 
perturbative effective superpotential to be of the form 



rpert _ w2 ^ I g kWt 1 m ikW* 1 
9i 2 (m T fk=1 )^T 



WlTf = WlF [ ^f^, - J / a k+i | (2.15) 



Since the theory depends on several masses g\, rrif k=1 , the function F can also depend 
on various ratios of them. From ( |2.15| ) we see that the number of W a in a diagram that 
contributes to W^Tf is constrained to be 

= 2 + J^ikr - 1) (2-16) 

r 

Here r labels the vertices of the Feynman diagram, and k r is the index k either for g k or 
for m* k at the vertex. We now write the diagram using 'tHooft's index loops notation. 
Every $ has a double line and every Q or Q has a single line. The number of index loops 
in the diagram is 

L = 2 + - ^2(k r - 1) - 2h - b (2.17) 

r 



where h is the number of handles in the surface and b is the number of boundaries (b is 
the number of Q or Q index loops). From ( |2.16| ) and ( |2.17| ) 

#W a = 2L-2 + 4h + 2b (2.18) 



As in [^6 1 , it is straightforward to extend this discussion to the case where the gauge group 



is broken, and the external light fields in the diagram are Wi l) . Here we find 

J2#W^ = 2L- 2 + 4/i + 26 (2.19) 

i 

(i) 

We also know that chiral operators have at most two factors of Wa in a trace and therefore 
there are at most two Wa^ on each index loop: 

J2#Wi i] <2L (2.20) 



Combining ( 2.19|) and ( |2.20| ) we learn that only three kinds of diagrams can contribute to 



the effective superpotential: 

1. h = b = 0. These diagrams have the topology of a sphere. One index loop has 
no operator insertion (leading to a factor of Ni) and every other index loop has an 
insertion of Sj leading to a result proportional to NiS^ -1 (S^ -1 stands for a product 
of L — 1 factors of Sj not necessarily all with the same value of j). 

2. h = b = 0. These diagrams also have the topology of a sphere. Two index loops have 

(i) 

an insertion of w a on each, and every other index loop has Sj leading to a result 
proportional to w a ^w^ a S^~ 2 . 

3. h = 0, b = 1. These diagrams have the topology of a disk. They have an insertion of 
Sj on each index loop leading to a result proportional to Sf. 

The sphere diagrams do not involve Q and Q. Therefore, their contribution is the 
same as in the theory without these fields; i.e. having only the adjoint field $. This 
theory has a shift symmetry associated with the decoupled diagonal U(l) field J2i w a 
This symmetry is best implemented by adding the auxiliary spinor coordinate ip a , 



and combining Si, Wa and Ni to a superfield 

Si(1>) =Si + f M i'» - V> V^z (2.21) 

Then the symmetry of shifting tfj a determines the sphere contribution in terms of a single 
function T§ ert of Si(t(j) 

' ' d 2 ^ r \SM)) (2.22) 
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The integral over ijj leads to the two sphere contributions we mentioned above. The function 
JFq 67 ^ is exactly the same as this function in the theory without the fundamental matter. 
It is given by 



2 ' *>J 

i¥=3 



Here we wrote explicitly the one loop result which depends on the UV cutoff Ao . The first 
term arises from integrating out the massive chiral superfield in the adjoint of U(Ni), and 
the second term is from the massive vector superfields which acquired mass in the Higgs 
mechanism. The higher order terms O(Sj) are independent of A . 

The disk diagrams lead to another function of Si whose contribution to the effective 
superpotential is not integrated over if; 

Tl ert =S> log (det ^) + <D(S*) (2.24) 

where we wrote explicitly the one loop terms which arise from the massive quarks Q and 
Q. Again, the higher order terms O(S'j) are independent of the UV cutoff A . jFf er * is not 
integrated over ip and it does not respect the shift symmetry of ip. This is consistent with 
the fact that in the theory with fundamental matter the over all U(l) superfield J2i w a^ 
does not decouple. 

The full effective superpotential at the scale fi where the gauge dynamics of the un- 
broken Y\i U(Ni) gauge fields is still weak is 

W e7r = ^E^l°g(^) 3 E / d2 ^fW + J **JF rt {SiU>))+lT rt {S i ) (2.25) 



The first term is the bare coupling. The second term includes the one loop running of the 
gauge couplings in Yl i SU(Ni); it depends only on the SU(Ni) fields Sj = Si + -^(wa^) 2 

J d^SfW=J2N tSl lo g ^ 3 (2.26) 
The dependence on the UV cutoff A from all the terms is 

(2iV-iV / )logAo^^ + logAo (S>« } ) ( 2 - 27 ) 



In accordance with the one loop beta function we choose the bare coupling 



r = -(2iV-iV / )logAo/A (2.28) 

with a finite scale A. This has the effect of changing Aq — > A everywhere except the term 
logA ■ This indicates that the over all U(l) superfield Yji w $ * s ^ ree i n ^ ne 

renormalized theory. We will continue to use the same notation as in ( |2.23|) - (|2.26|) but 
with Aq — > A. 

The strong IR dynamics is implemented by replacing (|2.26|) (after Aq — > A) with the 
Veneziano-Yankielowics superpotential P3j 

Wvy(s.) = J> ( 1o § ^m-+ N i)=lJ2j d ^ s " ( lo § ^ + 1) ( 2 - 29 ) 

Therefore, the final answer for the effective superpotential is 

W eff (S i ,w®,N i ) = J d 2 iPMSiW)+Fi{Si) = j d 2 ip (Jo(^))+V>Wi(W))) 

(2.30) 

with 
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(2.31) 

T ^pert 
J- 1 — J- 1 

Here = Si + ^jfiw^a) 2 are independent chiral superfields, and are independent field 
strengths of the massless vector superfields. 

Using the effective superpotential we can compute expectation values of operators by 
differentiating with respect to the coefficients in (|1.2p 



fc + 1 dg k dg k J dg k 

(M*- V> = ^ = — (2 ' 32) 

In the second equation we used the fact that the sphere diagrams do not involve the 
fundamental matter, and therefore the function JF is independent of rrij- k . The sphere 
contribution in the first equation is as in the theory without the fundamental matter. The 
disk contribution depends on mt k . Therefore (Tr $ fc+1 ) and (T(z)) depend on k in 
accord with the comment after (|2.10| ). 



Since we have already solved for these expectation values using the anomaly equations, 
we effectively found the full effective superpotential. The variables Si or more precisely Sj 
capture the information in the coefficients in the polynomial f(z) in (|2.5|) . Equivalently, 
the freedom in f(z) is determined by the values of the fields s,. The polynomials p a ( |2.10| ) 
are determined in terms of the variables w& . Finally, c(z) ( |2.10| ) is determined in terms 
of iVj. We remarked above that M^' (|2.14| ) are integrated out. This is closely related to 
the fact that the polynomials q(z) ( |2.8|) are determined (see the discussion after (|2.9|) ). An 
ambiguity in q(z) would have appeared as more fields in the effective superpotential. 



3. Matrix model 

We consider the matrix model with the "action" 

N 



A 

9 



Tr W(&) +QW f ($)Q 



(3.1) 



where $, Q and Q are iV x N, N x Nf and Nf x N dimensional matrices, and define its 
free energy T through 

exp f-^-Fj = J d^dQdQexp(-A) (3.2) 

Eventually, we will take the large iV limit with everything else held fixed. 
We will be interested in the resolvents 

R( z ) =t ( Tr 1 



N \ z-<5> 

M f M) = (Q7^Q f 

f \ 1 z - $ 

We perform the following three independent transformations 

5$ ' 



z - $ 



8Q 7 =Q~X 



1 

z — $ 

? 1 
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(3.3) 



SQ f = ^X f f ,Q f ' (3.4) 



to find 



N 
9 



z - $ 



z - $ 



Tr 



— Q— ^AQ = tr ATr ^ 

9 z - z — $ 

iV~~m(8)- ~ 1 

— QX^^Q = tr ATr ^ 

# z - $ 2 - $ 



(3.5) 



where we suppressed both the matrix indices and the flavor indices. Again we used the 
convention A f f B f ft , = (AB) f f '„, and tr denotes a trace over the flavor indices. Taking the 
expectation values of these equations and the large N limit they become 



W(z)R(z) = R{z 



(M(z)m(z)) f f ' 



R(z)5 f 



(m(z)M(z)) 



R(z)5 



f 
? 



(3.6) 



In the last two equations we used the fact that they are true for every A and A. 

The equations (|3.6|) are the loop equation |35[] of the matrix model. The first equation 
describes a closed string splitting to two closed strings. The other two equations describe 
an open string turning into a closed string. 

In solving the first equation for R(z) there is an ambiguity in a polynomial f(z) = 



-4 



W'(z)R(z) 



J + 



m(z)M(z) 



is 



However, in solving for M(z) the ambiguity in q(z) 

fixed by imposing that M(z) is not singular at the points where m(z) has a zero eigenvalue. 

The equations (|3l)|) are the same as the last three equations in the gauge theory ( |2~4l ) 
and the ambiguity in solving them (the value of f(z) and the way q(z) is determined) is 
also the same. Therefore we can identify 



R(z) 
M(z) 



:(R(Z)) 

<M(z)) 



It is interesting that while (Tr -^^) in the matrix model is not identified with (Tr 

z — $ 



(3.7) 
i \ 



z — $ 



z — $ 

in the gauge theory but with (Tr the matrix model object (Q—^Q) is identified 



with its natural gauge theory counterpart (Qjz^Q)- 

We now relate the effective superpotential in the gauge theory (|2.30| ) to the free energy 
of the matrix model T. We will need the sphere and the disk contributions to the free 
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energy: JF = lim^^^ JF, T\ = lim^^^ ^(JF — JF ). JF is computed in the theory without 
the fundamental matter. It was conjectured by Dijkgraaf and Vafa and proven in ||26|1 



(see also [p2[ ) that in this theory JF = JFq (recall that .Fo in the gauge theory ( |2.22|)(|2.30|) 



is the same as in the theory without the fundamental matter). We now turn to the disk 
amplitudes. Following ( IT?) ) 

(Qfb k - 1 Q f ) = {Qf> k ~ l Q f ) (3.8) 



They can be computed in the matrix model and in the gauge theory as derivatives with 

f ( 



respect to the coefficients in mt(z 



dm f fk dm f fk 



(3.9) 



where we used ( |2.32p and the fact that in the matrix model JF is independent of mt k . 



We learn that T\ = T\ + &T\ with bT\ independent of k . Using the special case with 
only mj k=1 = mSj nonzero and m — > oo we easily learn that bT\ = 0. 
We conclude that ^ 

(3.10) 

This establishes the equivalence of the gauge theory and the matrix model for this case. 
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